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Abstract 

From Landau-Alexander-McTague theory and Monte-Carlo simulation results we 
show that the modulated liquid obtained by subjecting a colloidal system to a peri- 
odic laser modulation has long range bond-orientational order and non-zero shear 
rigidity. From infinite field simulation results we show that in the modulated liquid 
phase, the translational order parameter correlation function decays to zero expo- 
nentially while the correlation function for the bond-orientational order saturates 
to a finite value at large distances. 
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1 Introduction 



Consider a 2-d charge stabilized colloidal system subject to a one-dimensionally 
modulated stationary laser field (obtained by superposing two beams). Chowd- 
hury, Ackerson and Clark [1] found that when the wave-vector of the laser 
modulation is tuned to be half the wave- vector q Q at which the structure factor 
of the colloidal liquid shows its first peak, above a certain laser field inten- 
sity, the system freezes into a 2-D triangular crystal. This phenomena of laser 
induced freezing (LIF) has been studied subsequently by experiments [2-4], 
simulations [5- 7] and density functional calculations [8-10]. For field strengths 
below the value at which the system undergoes freezing, the laser field in- 
duces a density modulation with wave-vector q in the liquid, and this phase 
has been called modulated liquid in the literature. In this paper we show that 
the modulated liquid has rather interesting properties. Specifically the in- 
duced translational order generates a field conjugate to the bond-orientational 
order parameter. Hence the "modulated liquid" has a non-zero value of bond- 
orientational order parameter, and consequently a finite rigidity (shear mod- 
ulus). We first present a qualitative picture in terms of a generalized Landau- 
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Alexander-McTague [11] theory and substantiate the picture with results from 
a detailed Monte-Carlo simulation. 



2 Mean-field treatment 



From general symmetry grounds the coarse-grained free energy functional in 
two dimensions in the presence of an external field V e , coupled to one of the 
six density Fourier modes (pg ) characterizing a 2-D crystal has the form [12]: 



F=~2V eP(?i + 1 -r T J2\PG\ 2 

G 

+w T ]T pgPg'Pg" + ME IpgI 2 ) 2 IpgI 4 + ••• 

G+G'+G"=0 G G 

+ ^r 6 \iJ 6 \ 2 + u 6 \iJ 6 \ 4 + ... 

+7E \Pg\ 2 IMG x - tG y f + r 6 (G x + iG y )% (1) 

G 

where the summations run over all the wave-vectors forming the first shell 
of the reciprocal lattice vector of the triangular lattice into which the colloid 
freezes and ipQ is the bond-orientational order parameter. 



Notations: 

Pi = P& = P-g, 

P2 = P G r 2 = Pg 3 = P-g 2 = P-G 3 



Fig. 1. Density Fourier modes for 2-D triangular crystal. 

In the absence of external field, for large positive values of and r 6 , the 
free energy is minimum when all order parameters are zero (ie. for the liquid 
phase). If, as a consequence of a change in temperature, screening length or 
some other parameter, tt decreases much faster than r 6 , then because of the 
cubic term in equation (1) one gets a first order freezing transition into a 
crystalline phase. On the other hand, if r$ decreases much faster than tt, 
one gets a continuous transition to an orientationally ordered hexatic phase, 
characterized by nonzero values of i/jq and zero values of pg. 

The external field induces a non-zero even in the liquid phase. The non- 
zero value of pg leads to an effective field conjugate to the bond-orientational 
order through the coupling 7 in eqn. 1. So ip§ is also turned on as the ex- 
ternal field is applied. Since the external field modulation fixes the directions 
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in which the order develops, we can treat the order parameters as real num- 
ber to get the equilibrium phase diagram (the phase in the order parameters 
will be important in determining the elastic coefficients but not in determin- 
ing the stable equilibrium phase). The experimental systems and the earlier 
simulations referred to above [1-7] corresponded to volume fraction and salt 
concentration values such that the zero field freezing transition was first or- 
der; hence they presumably correspond to tt < r§. Even if r$ is much smaller 
than r T , such that the freezing mechanism at zero external field is a two stage 
transition with the first (continuous) transition being to the hexatic phase, 
the external field will immediately destroy the intervening hexatic phase. 

0.3 




Fig. 2. Values of the translational and orientational order parameters as a function 
of external field in Landau- Alexander-McTague theory 



1.5, u'rp = 0.5 



Choosing the coefficients tt — r% — 0.30, wt = —1/3, ut = uq 
and 7 = —0.7 numerical minimization of eqn. 1 gives the order parameter pro- 
file as a function of external field shown in Fig. 2. i/jq and pg becomes nonzero 
as soon as V e is turned on, and at V e = 0.0239 all the order parameters jump 
simultaneously signifying a first order transition to a crystalline structure. 



One can define a rigidity modulus in analogy to the helicity modulus in su- 
perfluids [13]: 



Y = lim 



1 d 2 F 



q ^o q 2 d 2 q ' 



(2) 



where we consider an field coupled to the orientational order parameter as 
tp&iT) — > ip6{r)e lQ ' r ■ The free-energy cost for such a field will depend on the 
gradient term in (not included in eqn. 1). To leading order such a term 
would predict Y ~ V* as \q\ — > 0. 
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3 Monte-Carlo Simulations 



3.1 Simulation details 



Since the parameters in eqn. 1 are phenomenological and there is no direct 
way to fix them to characterize the experimental system, we have performed 
Monte-Carlo simulations to study the modulated liquid phase. We have con- 
sidered a 2-D system of charge stabilized colloidal particles (with diameter 
2R = 1.07/im) confined in a rectangular box of size ^a s L x a s L with peri- 
odic boundary conditions and subjected to an external potential of the form 
U(f) = —V e cos(q x)i with q = 27r/(^a s ), where a s is the mean inter-particle 
separation. The inter-particle interaction is modeled by the DLVO potential: 

tt (*.\ - (Ze) 2 f exp(KR) ^ exp(-Krij) 

Uij{r) ~ ^ { TT^R } K\ ( j 



Here Ze [Z = 7800) is the effective surface charge, e (=78) is the dielectric 
constant of the solvent and k is the inverse of the Debye screening length due to 
the small ions (counterions and impurity ions) in the solvent. The parameters 
are the same as in earlier simulation studies and similar to the experiments 

[!]■ 

To study the effect of very large external field (infinite field), in some of the 
simulations we had fixed the particles in parallel lines defined by the potential 
minima. Though the particles move freely only along the lines, they interact 
in full two dimensional space. The resulting simplification allows us to simu- 
late systems with L as large as 100, which allows us to compute correlation 
functions at large distances. 

The translational order parameters are defined by: 

Pc=<^E^>- (4) 



G\ refers to the direction parallel to the modulation wave- vector \ ^J 2a (1, 0)], 

while G 2 refers to the other independent wave-vector [ ^J 2a (§, ^ )] forming 
the first shell of wave-vectors of the triangular lattice. Value of as defined 
in eqn. (4) depends on the coordinate origin. So we measure the translational 



order parameter as pg =< -^y [J2i cosiG.fi)} 2 + siniG.fi)} 2 > [14]. This 
quantity is of order unity in crystalline state, while in liquid this goes to zero 
as ^t= for a system of N particles. 



4 



We define the bond-orientational order parameter as: 



1 1 

Jy m n 



>, 



(5) 



where B m ,n is the angle made by the line joining the position of particle m to 
the neighboring particle n, measured with respect to some fixed direction. z m 
is the number of neighbors corresponding to the Voronoi cell of particle m. 

To calculate order parameter correlation functions we define local order pa- 



rameters Pg(r^) = ^^i[e iG ' r " 1 + J2n eiG ' r "] and V^tjm) = j^J2n e(ii6m ' n , where 



the order parameters are defined at the position of a particle m, and n de- 
notes the particles forming the Voronoi cell of particle m. The order parameter 
correlation functions are defined as, 



To measure the helicity modulus Y, we have applied "anti-periodic" bound- 
ary condition along the x direction. In case of the standard periodic boundary 
condition, one repeats the simulation cell throughout space for calculating the 
inter-particle potential and to put back the particle inside the simulation box 
once it moves out in the course of the simulation. In case of the " anti-periodic" 
boundary condition, in the x direction for example, successive imaginary re- 
peat boxes are shifted by half the lattice spacing in the y-direction. Equiva- 
lently, while folding back a particle that exits from the simulation box in the 
±x direction, a displacement ±|a s y is applied. The change in energy between 
periodic and "anti-periodic" boundary conditions gives a measure of Y. 

3.2 results 

The specific heat in figure 3 (a) shows a peak as function of f3V e signifying 
a phase transition at f3V e = 0.25. In figure 3(b) we present order parameters 
for L = 20 and Ka s = 15.5 (where the system is liquid for zero external field) 
as a function of f3V e . The translational order parameters p\ and p 2 are non- 
zero even in the liquid phase because of the finite system size. Also they seem 
to grow continuously. But a finite-size scaling analysis shows that the results 
are consistent with first order transition scenario with small discontinuity in 
energy [7]. In addition to p±, ^ also shows non-zero value as soon as the 
external field is switched on as expected from the coarse-grained free energy. 

In figure 4, we plot p 2 and ^ for infinite field and 10 4 particle system (L = 




(6) 
(7) 
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Fig. 3. (a) Specific heat, and (b) Order parameters, obtained from Monte-Carlo 
simulations of a 400 particle system as a function of (3V e . na s was chosen to be 15.5 
for these simulations, such that the system is in liquid state in absence of external 
field 
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Fig. 4. and ^6 across the melting transition for infinite field and 10000 particles 

100) as a function of Ka s . While the translational order parameter shows a 
sharp fall at the melting transition (na s = 15.6), the bond-orientational order 
parameter remains large and finite. 



In figure 5 we plot the order parameter correlation functions in the crystalline 
region(5.a) and in the liquid region (5.b). In crystalline region Gt{d) decays 
as a power law, while Ge(y) saturates to a constant value. In the liquid region, 
while Gt(v) decays to zero exponentially at large distances, G&{y) still remain 
finite. 

In figure 6 we have shown the difference in energy in units of k^T per parti- 
cle for anti-periodic and periodic boundary conditions. We have considered a 
system with L = 10 in order to have a measurable enegy difference because 
of the large shear. The system freezes around (3V e = 0.2. From figure 6 we 
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Fig. 5. Translational and bond-orientational correlation functions along the external 
field minima (y axis) from simulation of 10000 particles and infinite field 
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Fig. 6. Energy difference between anti-periodic and periodic boundary conditions 

find that there is a non-zero elastic energy cost and hence a non-zero shear 
modulus even before the system freezes to crystalline structure. At f3V e = 0.2, 
there is a sharp change in AE, which then saturates. 



4 Conclusion 



We have shown that the modulated liquid phase obtained by inducing density 
modulation in a colloidal liquid by subjecting it to external laser modula- 
tions has properties intermediate between crystal and liquid; specifically par- 
tial translational order, finite bond-orientational order, and non-zero rigidity 
modulus. In two dimensional freezing, a large part of the loss of entropy upon 
freezing is due to the bond-orientational ordering. In LIF, the fact that the 
freezing transition occurs from a partially bond-orientationally ordered phase 
to a crystalline phase, helps one to understand why the transition becomes 
more weakly first order as the external field strength is increased. 
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